Elementary simulation of tethered Brownian motion

We describe a simple numerical simulation, suitable for an undergraduate project (or graduate problem set), of the Brownian motion of a particle in a Hooke-law potential well. Understanding this physical situation is a practical necessity in many experimental contexts, for instance in single molecule biophysics; and its simulation helps the student to appreciate the dynamical character of thermal equilibrium. We show that the simulation succeeds in capturing behavior seen in experimental data on tethered particle motion.Comment: Submitted to American Journal of Physic

DNA Looping Kinetics Analyzed Using Diffusive Hidden Markov Model

Tethered particle experiments use light microscopy to measure the position of a micrometer-sized bead tethered to a microscope slide via a ~micrometer length polymer, in order to infer the behavior of the invisible polymer. Currently, this method is used to measure rate constants of DNA loop formation and breakdown mediated by repressor protein that binds to the DNA. We report a new technique for measuring these rates using a modified hidden Markov analysis that directly incorporates the diffusive motion of the bead, which is an inherent complication of tethered particle motion because it occurs on a time scale between the sampling frequency and the looping time. We compare looping lifetimes found with our method, which are consistent over a range of sampling frequencies, to those obtained via the traditional threshold-crossing analysis, which vary depending on how the raw data are filtered in the time domain. Our method does not involve such filtering, and so can detect short-lived looping events and sudden changes in looping behavior.Comment: 3 page pdf including 3 figures corrections: 2nd page, 1st column, values of diffusion coefficient, spring constant and the decay time were typed incorrectly. No conlcusions were affecte

Tethered Particle Motion as a Diagnostic of DNA Tether Length

The tethered particle motion (TPM) technique involves an analysis of the Brownian motion of a bead tethered to a slide by a single DNA molecule. We describe an improved experimental protocol with which to form the tethers, an algorithm for analyzing bead motion visualized using differential interference contrast microscopy, and a physical model with which we have successfully simulated such DNA tethers. Both experiment and theory show that the statistics of the bead motion are quite different from those of a free semiflexible polymer. Our experimental data for chain extension versus tether length fit our model over a range of tether lengths from 109 to 3477 base pairs, using a value for the DNA persistence length that is consistent with those obtained under similar solution conditions by other methods. Moreover, we present the first experimental determination of the full probability distribution function of bead displacements and find excellent agreement with our theoretical prediction. Our results show that TPM is a useful tool for monitoring large conformational changes such as DNA looping

Calibration of Tethered Particle Motion Experiments

The Tethered Particle Motion (TPM) method has been used to observe and characterize a variety of protein-DNA interactions including DNA loping and transcription. TPM experiments exploit the Brownian motion of a DNA-tethered bead to probe biologically relevant conformational changes of the tether. In these experiments, a change in the extent of the bead’s random motion is used as a reporter of the underlying macromolecular dynamics and is often deemed sufficient for TPM analysis. However, a complete understanding of how the motion depends on the physical properties of the tethered particle complex would permit more quantitative and accurate evaluation of TPM data. For instance, such understanding can help extract details about a looped complex geometry (or multiple coexisting geometries) from TPM data. To better characterize the measurement capabilities of TPM experiments involving DNA tethers, we have carried out a detailed calibration of TPM magnitude as a function of DNA length and particle size. We also explore how experimental parameters such as acquisition time and exposure time affect the apparent motion of the tethered particle. We vary the DNA length from 200 bp to 2.6 kbp and consider particle diameters of 200, 490 and 970 nm. We also present a systematic comparison between measured particle excursions and theoretical expectations, which helps clarify both the experiments and models of DNA conformation

Diffusive hidden Markov model characterization of DNA looping dynamics in tethered particle experiments

In many biochemical processes, proteins bound to DNA at distant sites are brought into close proximity by loops in the underlying DNA. For example, the function of some gene-regulatory proteins depends on such DNA looping interactions. We present a new technique for characterizing the kinetics of loop formation in vitro, as observed using the tethered particle method, and apply it to experimental data on looping induced by lambda repressor. Our method uses a modified (diffusive) hidden Markov analysis that directly incorporates the Brownian motion of the observed tethered bead. We compare looping lifetimes found with our method (which we find are consistent over a range of sampling frequencies) to those obtained via the traditional threshold-crossing analysis (which can vary depending on how the raw data are filtered in the time domain). Our method does not involve any time filtering and can detect sudden changes in looping behavior. For example, we show how our method can identify transitions between long-lived, kinetically distinct states that would otherwise be difficult to discern

Spatial and topological organization of DNA chains induced by gene co-localization

Transcriptional activity has been shown to relate to the organization of chromosomes in the eukaryotic nucleus and in the bacterial nucleoid. In particular, highly transcribed genes, RNA polymerases and transcription factors gather into discrete spatial foci called transcription factories. However, the mechanisms underlying the formation of these foci and the resulting topological order of the chromosome remain to be elucidated. Here we consider a thermodynamic framework based on a worm-like chain model of chromosomes where sparse designated sites along the DNA are able to interact whenever they are spatially close-by. This is motivated by recurrent evidence that there exists physical interactions between genes that operate together. Three important results come out of this simple framework. First, the resulting formation of transcription foci can be viewed as a micro-phase separation of the interacting sites from the rest of the DNA. In this respect, a thermodynamic analysis suggests transcription factors to be appropriate candidates for mediating the physical interactions between genes. Next, numerical simulations of the polymer reveal a rich variety of phases that are associated with different topological orderings, each providing a way to increase the local concentrations of the interacting sites. Finally, the numerical results show that both one-dimensional clustering and periodic location of the binding sites along the DNA, which have been observed in several organisms, make the spatial co-localization of multiple families of genes particularly efficient.Comment: Figures and Supplementary Material freely available on http://dx.doi.org/10.1371/journal.pcbi.100067

Effect of promoter architecture on the cell-to-cell variability in gene expression

According to recent experimental evidence, the architecture of a promoter, defined as the number, strength and regulatory role of the operators that control the promoter, plays a major role in determining the level of cell-to-cell variability in gene expression. These quantitative experiments call for a corresponding modeling effort that addresses the question of how changes in promoter architecture affect noise in gene expression in a systematic rather than case-by-case fashion. In this article, we make such a systematic investigation, based on a simple microscopic model of gene regulation that incorporates stochastic effects. In particular, we show how operator strength and operator multiplicity affect this variability. We examine different modes of transcription factor binding to complex promoters (cooperative, independent, simultaneous) and how each of these affects the level of variability in transcription product from cell-to-cell. We propose that direct comparison between in vivo single-cell experiments and theoretical predictions for the moments of the probability distribution of mRNA number per cell can discriminate between different kinetic models of gene regulation.Comment: 35 pages, 6 figures, Submitte

Lac repressor mediated DNA looping: Monte Carlo simulation of constrained DNA molecules complemented with current experimental results

Tethered particle motion (TPM) experiments can be used to detect time-resolved loop formation in a single DNA molecule by measuring changes in the length of a DNA tether. Interpretation of such experiments is greatly aided by computer simulations of DNA looping which allow one to analyze the structure of the looped DNA and estimate DNA-protein binding constants specific for the loop formation process. We here present a new Monte Carlo scheme for accurate simulation of DNA configurations subject to geometric constraints and apply this method to Lac repressor mediated DNA looping, comparing the simulation results with new experimental data obtained by the TPM technique. Our simulations, taking into account the details of attachment of DNA ends and fluctuations of the looped subsegment of the DNA, reveal the origin of the double-peaked distribution of RMS values observed by TPM experiments by showing that the average RMS value for anti-parallel loop types is smaller than that of parallel loop types. The simulations also reveal that the looping probabilities for the anti-parallel loop types are significantly higher than those of the parallel loop types, even for loops of length 600 and 900 base pairs, and that the correct proportion between the heights of the peaks in the distribution can only be attained when loops with flexible Lac repressor conformation are taken into account. Comparison of the in silico and in vitro results yields estimates for the dissociation constants characterizing the binding affinity between O1 and Oid DNA operators and the dimeric arms of the Lac repressor. © 2014 Biton et al

Integration host factor alters LacI-induced DNA looping

The integration host factor protein of Escherichia coli, which sharply bends DNA at specific sites and non-specifically compacts the bacterial genome, can also alter looping of DNA in an artificial system based on the lactose repressor protein of E. coli. In single molecule experiments, we show that both specific bending and non-specific compaction alter LacI-mediated looping of DNA. Our results highlight the subtle regulatory roles that proteins, which confer structure upon DNA, might have in controlling DNA transcription and other processes in which the conformation of DNA determines the binding and activity of processive enzymes